Noname manuscript No. 

(will be inserted by the editor) 



Dinesh Singh 



An Analytic Perturbation Approach 
for Classical Spinning Particle 
Dynamics 



Received: date / Accepted: date 



Abstract A perturbation method to analytically describe the dynamics of a 
classical spinning particle, based on the Mathisson-Papapetrou-Dixon (MPD) 
equations of motion, is presented. By a power series expansion with respect 
to the particle's spin magnitude, it is shown how to obtain in general form an 
analytic representation of the particle's kinematic and dynamical degrees of 
freedom that is formally applicable to infinite order in the expansion. Within 
this formalism, it is possible to identify a classical analogue of radiative cor- 
rections to the particle's mass and spin due to spin-gravity interaction. The 
robustness of this approach is demonstrated by showing how to explicitly 
compute the first-order momentum and spin tensor components for arbitrary 
particle motion in a general space-time background. Potentially interesting 
applications based on this perturbation approach are outlined. 

Keywords classical spinning particles ■ spin-gravity interaction • perturba- 
tion approach 



1 Introduction 

For many years, there has been on- going research into better understanding 
the motion of extended bodies in strong gravitational fields, as described 
within the framework of general relativity. Since it is reasonable to expect 
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that such bodies wiU possess spin angular momentum during their formation, 
it is critical to properly account for all interactions involving the coupling 
of their spin to the curved space-time background. This understanding has 
relevance, for example, in identifying the motion of rapidly rotating neutron 
stars in orbit around supermassive black holes, a topic of particular interest 
for the space-based LISA gravitational wave observatory pj in observing low- 
frequency gravitational radiation emitted from these type of sources. 

The first recorded attempt to solve this problem was formulated by Mathisson 
[2] , who obtained an interaction term in the form of a direct spin coupling to 
the Riemann curvature tensor. Papapetrou [3] extended this initial idea by 
modelling the spinning particle as a matter field confined within a space-time 
world tube. Even more sophisticated approaches have been put forward by 
other people [4], primarily to deal with the higher-order multipole moment 
contributions to the spinning extended object's motion in curved space-time. 
In certain respects, it can be argued that Dixon [511S] produced the most 
complete model to date, which accounts for all multipole moment interac- 
tions to infinite order. For most practical purposes, however, it is sufficient 
to truncate the equations of motion to leading order in the spin, known as 
the "pole-dipole approximation" introduced by Mathisson and Papapetrou, 
on the condition that the the spinning object's dimensions are small com- 
pared to the local radius of curvature for the space-time background. Within 
this approximation scheme, these equations are commonly referred to as the 
Mathisson-Papapetrou-Dixon (MPD) equations. 

Besides a formal analysis of the MPD equations [IlIHllS] , there are partic- 
ular studies on spinning particle dynamics in a Kerr background [IO | ll H [T2 l 
113] , scattering interactions with gravitational waves |,14iil5j , and other appli- 
cations. In particular, the MPD equations lend themselves well to numerical 
analysis with applications involving gravitational wave generation in a Kerr 
background [I6ul7j . evidence of deterministic chaos within particle orbital 
dynamics p!8lfT9[|2Dll2T] . and particle motion in a Vaidya background with 
radially infalling radiation Recently, a first-order perturbative analysis 
of the MPD equations was developed by Chicone, Mashhoon, and Punsly 
(CMP) [23j, who applied their formalism to the study of rotating plasma 
clumps in astrophysical jets. As well, Mashhoon and Singh ,24j applied the 
CMP approximation to the orbital dynamics of spinning particles in a Kerr 
background, to compute leading order perturbations about a circular orbit 
and also explore the gravitomagnetic clock effect for spinning particles. 

This paper presents a perturbation analysis of the MPD equations follow- 
ing the initial attempt made by CMP, but now extended to arbitrary order 
in the perturbation expansion due to spin. For what is presented below, ge- 
ometric units of G = c = 1 are assumed and the metric has -1-2 signature. It 
begins with a brief outline of the MPD equations, followed by the perturba- 
tion approach adopted here. Some formal calculations of relevant kinematic 
and dynamical quantities based on this approach are then presented, leading 
to explicit calculation of the spinning particle's first-order momentum and 
spin tensor components for a particle with unspecified motion in a general 
space-time background. After discussion of some potential applications that 
may be relevant to this formalism, a brief conclusion follows. 
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2 Mathisson-Papapetrou-Dixon (MPD) Equations 

Equations of Motion: The starting point for describing spinning particle mo- 
tion in the "pole-dipole approximation" consists of the MPD equations of 
motion 

n 1 



dr 



- P'^ , (2) 
dr 

where P^^{t) is the spinning particle's linear momentum, S°'^{t) is the par- 
ticle's antisymmetric spin angular momentum tensor, u^{t) — da;''(T)/dT is 
the particle's centre-of-mass four-velocity, and R^vafi is the Riemann curva- 
ture tensor. While r will later become identified with the proper time of the 
particle along its centre-of mass worldline, within the context of the MPD 
equations it is strictly just a parametrization variable whose properties need 
to be specified separately. When considering more sophisticated models of 
spinning objects beyond the pole-dipole approximation, ^ and ^ will have 
extra terms of the form and T°'^ [711^ to denote additional forces and 
torques, respectively, based on higher-order multipole moments beyond the 
mass monopole and spin dipole moment. These extra terms require specifi- 
cation of the object's energy- momentum tensor T^'' [MSl llOll^l] . subject to 
covariant energy- momentum conservation T'^'^-i, = 0, which requires detailed 
knowledge of the particle's internal structure. However, for most practical 
purposes, the pole-dipole approximation is satisfactory. 

Supplementary Equations: By themselves, ([T]) and ([2]) are insufficient to com- 
pletely specify the motion of a spinning particle in curved space-time. To 
begin, the mass and spin parameters m and s are naturally identified in the 
form 



m 
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-P^. P"" , (3) 



«2 _ _ c qtii' 



(4) 



Furthermore, a supplementary spin condition needs to be specified in order to 
determine the particle's centre-of-mass trajectory. This is best accomplished, 
following Dixon's approach [51[B], by imposing an orthogonality condition 
between the particle's linear and spin angular momenta, in the form 

(5) 



While both ^ and ^ are formally functions of t, it can easily be shown [53] 
using the MPD equations and ^ that m and s are constants of the motion. 

With (O now specified, it is known [12] that the four-velocity m'' can be 
explicitly described in terms of P^ and S"^ , such that 



P ■ u 



(6) 
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where P ■ u is currently an undetermined quantity. A value for this scalar 
product needs to be chosen in order to determine the rate of the particle's 
internal clock with respect to r. It is important to emphasize from ([6]) that 
the linear momentum and four-velocity are not co-linear due to a non-trivial 
spin-curvature interaction. This has the effect of pulling the particle off a 
geodesic worldline, leading to interesting consequences for understanding the 
interplay between the particle's centre-of-mass motion and the dynamical 
response generated by the spin interaction with space-time curvature. 



3 Perturbation Approach to MPD Equations 

3.1 CMP Approximation 

As noted earlier, a first attempt in describing the MPD equations as a pertur- 
bation expansion was performed |231I24| , in which the underlying assumption 
is that — niu'^ = is small, where denotes the spin-curvature force. 
Furthermore, the spin magnitude is such that the M0ller radius p [23ll24l [25] 
is p = s/m <^ r, where r is the distance from the particle to the source gener- 
ating space-time curvature. This leads to the CMP approximation, expressed 
in the form 

D P'^ 1 

(7) 



dr 2 



0, (8) 



dr 

where the associated spin condition is 

5^. « . (9) 

It is important to recognize from ([5]) that, to first-order in s, the spin 
tensor is parallel transported in space-time. When comparing the CMP ap- 
proximation with a numerical integration of the MPD equations for circular 
motion around a Kerr black hole [24], it is shown that the kinematic prop- 
erties of the spinning particle agree very well in general. However, there is 
some loss of agreement when considering the radial component of the parti- 
cle's motion, since the CMP approximation does not reveal any modulation 
of its radial position compared to the corresponding plots generated by the 
MPD equations. This discrepancy appears when s/{mr) ^ 10~^ — 10~^ for 
r — 10 Af, where M is the mass of the Kerr black hole. It is most likely due 
to the lack of a more complicated spin interaction beyond what can be gen- 
erated by dS]). With this in mind, it seems appropriate to consider whether 
a more detailed approximation can account for the extra structure missing 
in the radial plots based on the CMP approximation alone, which is part of 
the motivation for this investigation. 



5 



3.2 Formalism 



Equations of Motion: The underlying assumption within this formalism is 
to define the particle's linear momentum and spin angular momentum in the 
form 



OO 

S'^'^{e)^eY,e^S^;^, (11) 

3=0 

where e is an expansion parameter corresponding to powers of s present in 
each order of the expansion, and P^j^ and S'|^'^ are the jth-order contributions 
of the linear momentum and spin angular momentum, respectively. That is, 
the zeroth-order quantities in e correspond to the dynamics of a spinless 
particle, while the first-order quantities in e lead to the CMP approximation 
described by ([7]) and ([8]). This is confirmed when substituting (fTOl) and (fTTj) 
into 

— = -- i?^a/3 S'^^ie) , (12) 

QT Z 

^^2^^^ =2£P[°(£)Z.^1(£), (13) 

where for ([T3|) the square brackets denote antisymmetrization of the indices 
defined by A^" B^^ = \ (y4" - Al^ B"), and an extra factor of e has to be 
added to maintain consistency with the definition (fTTj) . Furthermore, 

oo 

u^{e)^Y.^^ul^^, (14) 

3=0 

where the explicit expressions for u^^.^ are determined separately. When (llOp- 

([T4|l are combined together and grouped in terms of like orders of e, it follows 
that the jth-order expressions of the MFD equations are 

fc=0 



DS\ 



2^-2Vp[" .fl (16) 



k=0 



where, for P^^^ = mp u'^gs^ with 
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the zeroth- and first-order terms of (fTS]) and ^6]) in e are 



and 



dr 



- , (18) 



dr 



1 ^ ^ „^ 



. (20) 



This is consistent with the main results ([7]) and ([5]) of the CMP approxima- 
tion. 



Perturbation of the Mass and Spin Parameters: Besides the MPD equations, 
the supplementary equations (l3])-(l6]) serve a very important function in the 
development of this perturbation approach. This is especially true concerning 
the mass and spin magnitude parameters m and s, since they are dependent 
on P^ and S^'^, which are then expressible in terms of pH)) and pT|) . In 
fact, by invoking the language of quantum field theory, it becomes possible 
to identify the classical analogue of a bare mass toq defined according to (|17p . 
along with a bare spin sq, where 

such that the total mass and spin magnitudes exist as the sum of "radiative 
corrections" to mo and sq, due to the MPD equations in perturbative form. 
That is. 



'-ie)=e'sl\l + J2'' 



1 52 



(22) 
(23) 



where ruj and Sj are dimensionless jth-order corrections to mo and so, re- 
spectively, in the form 



-i = -iE^r'=^^'r.)> (24) 

fe=0 

1 ^ 

^ 72 X! '^Ai't '''' ^(k) ■ (25) 



^0 k=0 
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Explicit expressions of ([24]) and (|25|) required for subsequent calculations are 



and 



ml 



ml 



si 



,2 



(1) pM 



(0) 



O p(2) _L p(l) pM 



p(3) pM , p(2) 



0(1) CM" 



O c(2) QiJ-f I c(l) cMi' 
^ '^/ji^ "^(0) ^ '-^A"' (1) 



(26) 
(27) 
(28) 

(29) 
(30) 



Since m and s are already known to be constant for an exact treatment of 
the MPD equations, it follows that fhj and Sj must also be constant for each 
order of the perturbation expansion in e. While it is not obvious that this 
property emerges from ((24|) and ((25)) . this can be verified to at least the orders 
of expansion considered in this paper, upon evaluating the perturbative form 
of u^'. 



3.3 Kinematic and Dynamical Quantities 



Four- Velocity: As noted earlier, supplementary equations are needed to com- 
pletely specify the motion of a spinning particle according to the MPD equa- 
tions. This includes the four-velocity described by ^ once the undetermined 
scalar product P • u is specified. It turns out that by setting 



it follows that 
where 



m{£) 



P ■ u = — m(e), 

1 S^''{e)R,^^0 P-^je) 5°^(£)' 



1 



4 m2(e) 



m2(e)Z\(£) 



R^..o.0 S'^'^is) S'^^ie 



(31) 



(32) 



(33) 



With this particular choice of parametrization constraint, it is straightfor- 
ward to see that ([32|) . along with the spin condition constraint equation ([5|, 
leads to 



u^(e)u'^(e) = -1 + 



1 



■R^{e)R''{e) = -1 + 0{e^) , (34) 
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where 



(35) 

At least to third-order in e, this expression for ([34]) assuming (|3ip justifies 
the labeUing of t as proper time for parameterizing the spinning particle's 
motion along its centre-of-mass worldline. 

Though somewhat tedious to calculate, it is a straightforward exercise to 
evaluate the four- velocity u'^ perturbatively, based on p2|) and using (fTO)) . 
(fTTj) . and ([22l) . Therefore, it follows that the explicit expression for the spin- 
ning particle's four-velocity in general form is 



MO) 

mo 



1 



mo 



1 



— ( P{i) 2 ^(0) 



mo 



^(^2) - 2 " 2 



'1 I ^(0) 



mo 



^(3) - 2 '^i ^(2) " 2 



rn^m^-H-mnP(';) 



2m5 



Qa/3 - 2 Q/^J^ p7 cct/^ 



.n=0 



fc=0 



Given 
form 



(36) 



, along with the spin condition constraint equation ([5]) in the 



p(0) cM"^ _ _ p(fe) CM" 



fc) 



J> 1 



fc=i 



for the (j-l-l)th-order contribution in e, with 

"^(O) — " ' 

p(0) oMi' I p(l) cM'^ _ n 

'5(1) + -Sj-Q^ - U , 



(37) 

(38) 
(39) 
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for the first- and second-order perturbations in e, respectively, it is straight- 
forward to verify that p6|) satisfies (pi]) to third order in e. It follows that 
the perturbation of the spinning particle's worldline away from a geodesic is 
attainable by integration of ([55)1 with respect to r. 

Constancy of the Mass and Spin Parameters: With (|36|) evaluated, it is pos- 
sible to verify that 

^ - ^ ^ (40) 



dr dr 

for at least (|26|) - ((30)) . given (fT5| and p6| . In particular, from evaluating 
Dm^/dr directly according to p^ . it is found that 



= R,.o.p V V ur,_,_,, S^i'. , (41) 




which is identically satisfied for the required expressions of u'^j^ given by 

(|36p . Furthermore, it can be shown explicitly that the first- and second-order 
perturbations of the spin tensor satisfy 



dr 



0, (42) 



(£1 — _ p^^' cH'^ n p7 

- ^3 ^(0) -^(0) ^'^1^13 -f(o) '3(0) ■ 

M0ller Radius: An important quantity to evaluate via this perturbation ap- 
proach is the M0ller radius p — s/m^ which has significance in determining, 
for example, the strength of the spin-curvature force when applied to circular 
orbits of spinning particles in black hole space-times [221124] . Getting a bet- 
ter detailed sense of how the M0ller radius appears due to this perturbation 
approach may become useful for understanding the precise conditions for a 
transition into chaotic dynamics, as suggested in previous work [T8 l ll9 1 [2 m 
[21] . Given (|26p - p0)) . it is a straightforward calculation to show that 



m(e) Too 



2 ^ ^ 



3mt) 



+ 0(£4) ,(44) 



where non-trivial deviations from po = sq /toq due to geodesic motion emerge 
at second order in e. Since the "radiative corrections" to too and So described 
by ([M]l and (P5)l have the effect of increasing both the particle's mass and 
spin simultaneously, the M0ller radius appears to remain in the vicinity of pq. 
At present, however, it is impossible to precisely determine the nature of the 
shift from po to p without first examining (I44p with respect to a particular 
space-time background. 



10 



3.4 First-Order Linear Momentum and Spin Components 

Local Fermi Co-ordinate System: To fully appreciate the value of this per- 
turbation approach for describing the MPD equations, it is useful to compute 
the first-order linear momentum and spin angular momentum components in 
general form. Evaluation of P^'^j is particularly straightforward when formu- 
lating the problem in terms of a local Fermi co-ordinate system [24] and an 
orthonormal tetrad frame A'^q,. This leads to Pl^^^ — X^^a ^(1)7 where a denote 
indices for the Fermi co-ordinates x" defined on a locally flat tangent space 
in the neighbourhood of the spinning particle, and P"^-^ is the corresponding 
expression for the linear momentum on the tangent space. The tetrad frame, 
with A^Q = ''^('o)' satisfies the orthonormality condition 

V&$ = A^a A"^^ , (45) 
and satisfies the parallel transport law 

^^^^ (46) 



dr 



with respect to the general space-time co-ordinates x'^. The Riemann curva- 
ture tensor projected onto the locally flat tangent space then satisfies 

^^aM'S ^ ^t^'^P'^ ^''7 ■^'^5 ■ i^'^) 

In addition, the first-order spin condition with ((^ requires that S'^^^^ = 
A^j X^j S^ff^-^ to preserve orthogonality with P|,p^ in terms of the tetrad frame 

m. 



First- Order Linear Momentum: Given the tetrad formalism presented here, 
evaluation of the first-order linear momentum components becomes very easy 

to implement. From HI]) and ([46]), it follows that DP^'^^/dr = A'^s (dP^^^/dr^ , 

which results in 

^ = -2 ^%-^(0) (48) 
to be solved. This leads to the general expression 



Pii)--l>^'^ I (^^%-^|))dr. (49) 
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First-Order Spin Angular Momentum: Determining the first-order spin ten- 
sor components, in contrast to that of the linear momentum, is not so 
straightforward. This is because a similar line of reasoning when applied 
to leads to the expression 

-jiil = A-A--iil = 0, (50) 

in which the dS^f^/dr defined on the local tangent space are still undeter- 
mined. Specifically, it is not necessarily true that dS'"j'^/dr = 0, since the six 
equations in (jSOp may not all be linearly independent of each other. Further- 
more, even if this requirement were to be satisfied, the resulting constants of 

integration S"f^ cannot be further specified without some physical justifica- 
tion. Therefore, a different approach is required. 

Fortunately, there is a means available to solve for in general form, 
based on making full use of the spin condition constraint equation ([5]). Given 
, there are four equations of the form 

A^S^^'^-B-' = 0, (51) 

where 

A, ^ Plf^ , (52) 

B-^ ^ ^P^'^ S^,]. (53) 

As well, contraction of ([39]) with leads to two constraint equations of 
the form 

c^,5[;'; = o, (54) 

E^,, S^;'^ = , (55) 



where 

C,^^pI"^pI'\ (56) 
E,,^Pl"jDPl^^/dr, (57) 

given (UHl) and (l42|l . 

With (|5T1) . ((54| . and ((55| . there are six linear equations in S^^-^ and six 
unknowns for the spin tensor components, so it should be possible to solve 
for a unique solution of Sf^^^ algebraically. However, there is a complication 
in that one of the equations in (j5ip is linearly dependent on the other three, 
leading to the first-order spin condition constraint equation (|38p . Therefore, 
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the linear system of equations is formally underdetermined in one variable, 
leading to the solution 

\Hi2 Jl3 ~ J12 Hi3 I Aq \Hi2 Ji3 — J12 Hi3 I 



_ f cl3 _ ( o23 _ '^Oj 



^02 cl2 , f r,23 , 



^01 _ /■ ^2A r.12 , ( C.13 



expressed in terms of Sj^-^ , and 

=Cy+2C7oM,]Mo, (63) 

Jy =i;y+2£;o[,A,]Mo. (64) 
The expressions ((58l) - ([62|) can then be put into the form 

5^;'^ = M'"' 5g + Nf"" , (65) 
where M^''(r) and N'^''{t) are antisymmetric in their indices, and 

Af23 ^ 1 ^ iv23 = . (66) 

There is, however, one remaining equation to consider, namely (j42p which 
describes the parallel transport of S^^-^ . Since this is true for every component 

of the first-order spin tensor, it must certainly apply for the case of S'^^y 
Therefore, given (|42)) and ((65| , it follows that the first-order linear differential 
equation 

-^— + P{r)Sf!^{r)^Q{r) (67) 
must be satisfied, where 

Pir) = -2uf,^rj^pM'^^, (68) 
Q(r) = 2ufo)^i'^iV3l^ (69) 
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in terms of the metric connection -Tj^. Since the solution to (1671) is known 
exactly [26] in terms of a suitably chosen integrating factor /i(T), it follows 
that 

^a)M = ^ (^%(r')0(r')dr') , (70) 

where 

m(t) =cxp^ P(r')dT'^ , (71) 

and S'I^^^^(O) ~ 0. Given the generally non-trivial dependence of r on the in- 
tegrands in ((70|) and fTT]) . it is most likely that S'^^.j'^ will have to be solved 
numerically. While presented in a compact form, it is rather surprising that 
the solution ([70]) possesses a considerable level of complexity. It will be in- 
teresting to analyze the structure of S^^^ when calculated for particle motion 
in a simple but non-trivial space-time background, such as the Schwarzschild 
metric. 

4 Potential Applications 

It seems clear that this perturbation approach to the MPD equations leads 
to a robust formalism that lends itself well to various applications. For ex- 
ample, following upon previous work [24], it is possible to investigate the 
dynamics of a spinning particle in circular orbit around a Schwarzschild or 
Kerr black hole in vacuum [27] , or in the presence of radiation, as described 
by the Vaidya or Kerr-Vaidya metrics [28 . As mentioned earlier, having an 
analytic perturbative representation of the MPD equations may become very 
useful for identifying the conditions for a transition from stable motion to 
a chaotic form. In addition, it may be possible to study the spin-spin inter- 
action between two or more spinning particles, and determine the possible 
impact on tidal acceleration effects experienced by a reference particle within 
this configuration. There are likely many other applications to follow from 
future considerations of this formalism. 

5 Conclusion 

This paper has displayed a systematic expression of the Mathisson-Papapetrou- 
Dixon (MPD) equations in perturbative form that can be formally extended 
to infinite order in the perturbation expansion parameter e, corresponding to 
the order of the particle's spin magnitude s. If the spinning extended object 
in motion is treated as a test particle with a M0ller radius p = s/m <^ r, 
where r is the radius of curvature defined by the gravitational source, then 
the relevant expressions for the perturbation will converge rapidly, and only 
the first- or second-order quantities in e are likely required for most practical 
calculations. It may be useful to consider how to extend this formalism for 
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a spinning object that does not satisfy the "pole-dipole approximation," the 
results of which may perhaps become a suitable barometer for comparing 
competing models of extended objects in curved space-time. This may be a 
topic worthy of a future investigation. 
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